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The Twente turbulent Taylor-Couette (T 3 C) facility: Strongly turbulent (multiphase) 
flow between two independently rotating cylinders 
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A new turbulent Taylor-Couette system consisting of two independently rotating cylinders has been 
constructed. The gap between the cylinders has a height of 0.927 m, an inner radius of 0.200 m and a 
variable outer radius (from 0.279 to 0.220 m). The maximum angular rotation rates of the inner and 
outer cylinder are 20 Hz and 10 Hz, respectively, resulting in Reynolds numbers up to 3.4 x 10 6 with 
water as working fluid. With this Taylor-Couette system, the parameter space (Re*, Re , n) extends 
to (2.0 x 10 6 , ±1.4 x 10 6 , 0.716 — 0.909). The system is equipped with bubble injectors, temperature 
control, skin-friction drag sensors, and several local sensors for studying turbulent single-phase and 
two-phase flows. Inner-cylinder load cells detect skin-friction drag via torque measurements. The 
clear acrylic outer cylinder allows the dynamics of the liquid flow and the dispersed phase (bubbles, 
particles, fibers etc.) inside the gap to be investigated with specialized local sensors and nonintrusive 
optical imaging techniques. The system allows study of both Taylor-Couette flow in a high-Reynolds- 
number regime, and the mechanisms behind skin-friction drag alterations due to bubble injection, 
polymer injection, and surface hydrophobicity and roughness. 
[doi |10.1063/l 3548924] 



I. INTRODUCTION 



A. Taylor-Couette Flow 



Taylor-Couette (TC) flow is one of the paradigmatical 
systems in hydrodynamics. It consists of a fluid confined 
in the gap between two concentric rotating cylinders. TC 
flow has long been known to have a similarity PHI] to 
Rayleigh-Benard (RB) flow, which is driven by a tem- 
perature difference between a bottom and a top plate in 
the gravitational field of the earth [J, |5j . Both of these 
paradigmatical hydrodynamic systems in fluid dynamics 
have been widely used for studying the primary insta- 
bility, pattern formation, and transitions between lami- 
nar flow and turbulence @. Both TC and RB flows are 
closed systems, i.e., there are well-defined global energy 
balances between input and dissipation. The amount of 
power injected into the flow is directly linked to the global 
fluxes, i.e. angular velocity transport from the inner to 
the outer cylinder for the TC case, and heat transport 
from the hot bottom to the cold top plate for RB case. 
To obtain these fluxes, one only has to measure the cor- 
responding global quantity, namely the torque required 
to keep the inner cylinder rotating at constant angular 
velocity for the TC case, and the heat flux through the 
plates required to keep them at constant temperature for 
the RB case. In both cases the total energy dissipation 
rate follows from the global energy balances Q. From 
an experimental point of view, both systems can be built 
with high precision, thanks to the simple geometry and 
the high symmetry. 

For RB flow, pattern formation and flow instabilities 
have been studied intensively over the last century in low- 
Rayleigh (Ra) numbers, see e.g. review Q- With increas- 
ing Rayleigh numbers, RB flow undergoes various tran- 



sitions and finally becomes turbulent. In the past twenty 
years, the investigation on RB flow has been extended to 
the high-Rayleigh-number regime, which is well beyond 
onset of turbulence. To vary the controlled parameters 
experimentally, RB apparatuses with different aspect ra- 
tios and sizes have been built in many research groups in 
the past twenty years [8l-[l6j . Direct numerical simulation 
(DNS) of three-dimensional RB flow allows for quantita- 
tive comparison with experimental data up to Ra = 10 11 
17], which is well beyond the onset of turbulence 
The dependence of global and local properties on control 
parameters, such as Rayleigh number, Prandtl number, 
and aspect ratio, has been well-explored. The RB sys- 
tem has shown surprisingly rich phenomena in turbulent 
states (see, e.g., the review articles in Refs. [4( and [5|), 
and it is still receiving tremendous attention (see Ref. 
.18]). 

With respect to flow instabilities, flow transitions, and 
pattern formation, TC flow is equally well-explored as RB 
flow. Indeed, TC flow also displays a surprisingly large 
variety of flow states just beyond the onset of instabili- 
ties [19l - l21j . The control parameters for TC flow are the 
inner cylinder Reynolds number Rej, the outer cylinder 
Reynolds number Re G , and the radius ratio of the inner to 
outer cylinders rj = r%/r . Similar to RB flow, TC flow 
undergoes a series of transitions from circular Couette 
flow to chaos and turbulence with increasing Reynolds 
number (22[. The flow-state dependence on the rotation 
frequencies of the inner and outer cylinder in TC flow 
at low Reynolds numbers has been theoretically, numeri- 
call y a nd experimentally well-studied in the last century 
[23l428j . This is in marked contrast to the turbulent case, 
for which only very few studies exist, which we will now 
discuss. 

Direct numerical simulation of TC flow is still limited 
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to Reynolds numbers up to 10 4 , which is still far from 
fully-developed turbulence 0, [3(| ■ Experimentally, few 
TC systems are able to operate at high Reynolds num- 
bers (Rei > 10 4 ), which is well beyond the onset of 
chaos. Smith and Townsend (3lL l32j performed velocity- 
fluctuation measurements with a hot-film probe in turbu- 
lent TC flow at Rei ~ 10 4 , when only the inner cylinder 
was rotating. Systematic local velocity measurements on 
double rotating systems at high Reynolds numbers can 
only be found in Ref. [33| from 1933. However, the mea- 
surements [33j were performed with pitot tubes, which 
are an intrusive experimental technique for closed sys- 
tems. Recently, Ravelet et al. built a Taylor- Couette 
system with independently rotating cylinders, capable of 
Reynolds numbers up to 10 5 , and they performed flow 
structure measurements in the counter-rotation region 
using particle image velocimetry. 

The most recent turbulent TC apparatus for high 
Reynolds numbers (Rei ~ 10 6 ) was constructed in Texas 
by Lathrop, Fineberg and Swinney in 1992 [3^, [36| . 
This turbulent TC setup had a stationary outer cylin- 
der and was able to reach a Reynolds number of Re = 
1.2 x 10 6 . Later, the system was moved to Lathrop's 
group in Maryland. This setup will be referred to as the 
Texas-Maryland TC, or T-M TC in short. The appara- 
tus was successfully used to study the global torque, local 
shear stress, and liquid velocity fluctuations in turbulent 
states [35^40]. 

The control parameter Rei was extended to Rei ~ 10 6 
by the T-M TC; however, the roles of the parameters Re Q 
and 77 in the turbulent regime still have not been studied. 
Flow features inside the TC gap are highly sensitive to 
the relative rotation of the cylinders. The transition path 
from laminar flow to turbulence also strongly depends on 
the rotation of both the inner and the outer cylinders. A 
system with a rotatable outer cylinder clearly can offer 
much more information to better understand turbulent 
TC flow, and this is the aim of our setup. On the theo- 
retical side, Eckhardt, Grossmann, and Lohse [3j extend 
the unifying theory for scaling in thermal convection (4ll — 
l44j ] from RB to TC flow, based on the analogy between 
RB and TC flows. The gap ratio is one of the control pa- 
rameters in TC flow, and it corresponds to the Prandtl 
number in Rayleigh-Bcnard flow [3j. They define the "ge- 
ometrical quasi-Prandtl number" as {[(1 + 77) / '2] / y/rj} 4 ■ 
In TC flow the Prandtl number characterizes the geome- 
try, instead of the material properties of the liquid as in 
RB flow Q . From considerations of bounds on solutions 
to the Navier-Stokes equation, Busse [45[ calculated an 
upper bound of the angular velocity profile in the bulk 
of the gap for infinite Reynolds number. The predic- 
tion suggests a strong radius ratio dependence. It is of 
great interest to study the role of the "geometrical quasi- 
Prandtl number" in TC flow. This can only be done in 
a TC system with variable gap ratios. 

The T-M TC system was designed 20 years ago and 
unavoidably exhibits the limitations of its time. To ex- 
tend the parameter space from only (Rei, 0, fixed if) to 



(Re^ Re D , variable 77), we built a new turbulent TC sys- 
tem with independently rotating cylinders and variable 
radius ratio. 



B. Bubbly drag reduction 

Another motivation for building a new TC system is 
the increasing interest in two-phase flows, both from a 
fundamental and an applied point of view. For exam- 
ple, it has been suggested that injecting bubbles under 
a ship's hull will lower the skin-friction drag and thus 
reduce the fuel consumption; for a recent review on the 
subject we refer to Ref. 46]. In laboratory experiments 
skin-friction drag reductions (DR) by bubble injection up 
to 20% and beyond have been reported (13, E^. However, 
when supplying an actual-scale ship with bubble genera- 
tors, the drag reduction drops down to a few percent [49| . 
not taking into account the power needed to generate the 
air bubbles. A solid understanding of the drag reduction 
mechanisms occurring in bubbly flows is still missing. 

The conventional systems for studying bubbly DR are 
channel flows [5(|, flat plates (UH^I, and cavitation flows 
(46| . In these setups it is usually very difficult to control 
the power input into the flow, and to keep this energy 
contained inside the flow. In 2005, the T-M TC was out- 
fitted with bubble injectors in order to examine bubbly 
DR and the effect of surface roughness [H, [HJ. The 
strong point of a Taylor- Couette system, with respect 
to DR, is its well defined energy balance. It has been 
proved that the turbulent TC system is an ideal system 
for studying turbulent DR by means of bubble injection 
[H-[5|. 

The previous bubbly DR measurements in TC flow 
were based only on the global torque, which is not suffi- 
cient to understand the mechanism of bubbly DR. Vari- 
ous fundamental issues are still unknown. How do bub- 
bles modify the liquid flow? How do bubbles move inside 
the gap? How do bubbles orient and cluster? What is 
the effect of the bubble size? These issues cannot be 



addressed based on the present TC system. Ref. [53 1 
found that significant DR only appears at Reynolds num- 
bers larger than 5 x 10 5 for bubbles with a radius ~ 1 
mm. The maximum Reynolds number of the T-M TC is 
around 10 6 , which is just above that Reynolds number. 
A system capable of larger Reynolds numbers is therefore 
favored for study of this hitherto-unexplored parameter 
regime. The influence of coherent structures on bubbly 
DR can be systematically probed with a TC system when 
it has independently rotating cylinders. 



C. Twente turbulent Taylor-Couette 

Using the design of the T-M TC system as a start- 
ing point, we now present a new TC system with in- 
dependently rotating cylinders, and equipped with bub- 
ble injectors, dubbed "Twente turbulent Taylor-Couette" 
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FIG. 1: A photograph of the T C system. The height of the cylinder is about 1 m. 



(T 3 C). We list the main features of the T 3 C facility (the 
material parameters are given for water at 21°C as the 
working fluid): 

• The inner and outer cylinder rotate independently. 
The maximum rotation frequencies for the inner 
and outer cylinder are /j = 20 Hz and f = 10 Hz, 
respectively. 

• The aspect ratio and radius ratio are variable. 

• The maximum Reynolds number for the counter- 
rotating case at the radius ratio of rj = 0.716 is 3.4 
x 10 6 . The Reynolds number for double rotating 
cylinders is defined as Re — ($7^ — f2 r )(r — rj)/V, 
where fli — 2ir fa and tt = 27r/ are the angular 
velocities of the inner and outer cylinder, and v is 
the kinematic viscosity of water at the operation 
temperature. 

• The maximum Reynolds numbers at the radius ra- 
tio of T] = 0.716 are Re^ = CliTi(r — r^jv = 
2.0 x 10 6 for inner cylinder rotation and Re c = 
^0^0(^0 — Ti)/v = 1.4 x 10 6 for outer cylinder rota- 
tion. 

• Bubble injectors are incorporated for injecting bub- 
bles in a range of diameters (100 fim to 5 mm), 
depending on the shear strength inside the flow. 



• The outer cylinder and parts of the end plates 
are optically transparent, allowing for optical mea- 
surement techniques like laser Doppler anemome- 
try (LDA), particle image velocimetry (PIV) and 
particle tracking velocimetry (PTV). 

• Temperature stability and rotation rate are pre- 
cisely controlled. 

• Local sensors (local shear stress, temperature, 
phase-sensitive constant temperature anemometry 
(CTA), etc.) are built in or are mountablc. 

Initial torque data from the T 3 C and an upgraded T- 
M experiment recently appeared in Refs. [56l. 157]. This 
paper serves to detail the T 3 C apparatus, capabilities 
and measurement techniques. 



II. SYSTEM DESCRIPTION 

Figure [T] is a photograph of the T 3 C mounted in the 
frame. The details of the system will be described in 
Sees. IHAWTI 
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A. Geometry and materials 

As shown in Fig. [2j the system contains two indepen- 
dently rotating cylinders of radii and r Q . The working 
liquid is confined in the gap between the two cylinders of 
width d — r a — r^. The height of the gap confined by the 
top and bottom plate is L. By design, one set of radius 
ratios (77 = T jvi) and aspect ratios r = L/(r Q — r-j) of the 
T 3 C nearly match those of the T-M TC in order to allow 
for a comparison of the results. To increase the system 
capacity for high Reynolds numbers, the present T 3 C sys- 
tem has twice the volume compared to that of the T-M 
TC. The maximum Reynolds number of the T 3 C system 
is 3.4 x 10 6 when the two cylinders at a radius ratio of 
rj = 0.716 are counter-rotating with water at 21 °C as the 
working fluid. Table U lists the geometric parameters. 




FIG. 2: Schematic Taylor-Couette setup consisting of two in- 
dependently rotating coaxial cylinders with angular rotation 
rates Qi and fi . The gap between the cylinders is filled with 
a fluid. 



TABLE I: Geometric parameters of both TC setups, with 
inner radius r», outer radius r D , gap width d, radius ratio rj, 
aspect ratio F and gap volume V gap - The outer radius of the 
T 3 C gap can be varied resulting in different aspect and radius 
ratios. The value of 0.220 in the brackets refers to the case of 
a sleeve around the inner cylinder. 





T-M TC 






C 




L (m) 


0.695 


0.927 








imid (m) 


0.406 


0.536 








n (m) 


0.160 


0.200 


(0.220) 






r a (m) 


0.221 


0.279 


0.260 


0.240 


0.220 


d — r a ~ n (m) 


0.061 


0.079 


0.060 


0.040 


0.020 


V = ri/r a 


0.725 


0.716 


0.769 


0.833 


0.909 


r = L/d 


11.43 


11.68 


15.45 


23.18 


46.35 


Vgap (m 3 ) 


0.051 


0.111 


0.080 


0.051 


0.024 



As shown in Fig. [2] the inner cylinder (IC) consists of 
three separate sections ICbot, IC m id and ICt op , each able 
to sense the torque by means of load cell deformation 
embedded inside the arms connecting the IC sections to 
the IC drive shaft. End effects induced by the bottom 
and top of the TC tank are significantly reduced when 
focusing only on the IC m id section. The gap between 
neighbouring sections is 2 mm. The material of the IC 
sections is stainless steel (grade 316) with a machined 
cylindricity (radial deviations) of better than 0.02 mm. 

The outer cylinder (OC) is cast from clear acrylic, pro- 
viding full optical access to the flow between the cylin- 
ders. The OC is machined to within tolerances by Blan- 
son Ltd. (Leicester, UK) and the final machining was per- 
formed by Hemabo (Hengelo, Netherlands), which con- 
sisted of drilling holes in the OC for sensors, and remov- 
ing the stresses in the acrylic by temperature treatment. 
The thickness of the OC is 25.4 mm. The bottom and 
top of the TC tank are connected by the OC and rotate 
as one piece, embedding the IC completely with the IC 
drive axle protruding through the top and bottom plates 
by means of mechanical seals. 

The frame itself, as shown in blue in Fig. [TJ is sup- 
ported by adjustable air springs underneath each of its 
four support feet; they lift the frame fully off the ground. 
In combination with an inclination sensor, the frame can 
automatically level itself to ensure vertical alignment of 
the IC and OC with respect to the gravity. 



B. Varying gap width and IC surface properties 

Additional clear acrylic cylinders are available for this 
new T 3 C system, and can be fitted between the orig- 
inal OC and IC. These "filler" cylinders will rotate to- 
gether with the original OC and provide a way to decrease 
the outer radius of the gap. The nominal gap width of 
0.079 m can thus be varied to gap widths of 0.060, 0.040 
and 0.020 m, resulting in the aspect and radius ratios as 
shown in Table [I] 

Instead of a smooth stainless steel IC surface, other IC 
surfaces with alternative chemical and/or surface struc- 
ture properties can be employed, preferably in a perfectly 
reversible way. Thus we employed a set of cylindrical 
stainless steel sleeves, installed around the existing IC 
sections, leaving the original surface unaltered. These 
sleeves each clamp onto the IC sections by means of a pair 
of polyoxymethylene clamping rings between the sleeve 
and the IC. The inner radius of the gap is hence increased 
by 0.020 m. For good comparison at least two sets of 
sleeves need to be available; one set with a bare smooth 
stainless steel surface to determine the effect of a chang- 
ing gap width, and one set with the altered surface prop- 
erties. The sleeves themselves can be easily transported 
for surface-altering treatment. 
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FIG. 3: A sketch of the signal and control system of the T 3 C system. The abbreviations in the sketch: O.C. (outer cylinder), 
I.C. (inner cylinder), DAQ (data acquisition), PLC (programmable logic controllers), and PC (personal computer). 



C. Wiring and system control 

A sketch of the control system is shown in Fig. [3] All 
sensors embedded inside the IC are wired through the 
hollow IC drive axle, and exit to a custom-made slip 
ring from Fabricast (El Monte, USA) on top, with the 
exception of the shear stress CTA (Constant Tempera- 
ture Ancmometry) sensors described in Sec. Ill II The 
slip ring consists of silver-coated electrical contact chan- 
nels with 4 silver-graphite brushes per channel. The 4 
brushes per channel ensure an uninterrupted signal trans- 
fer during rotation, and increase the signal-to-noise ra- 
tio. The most important signals being transferred are 
the signal-conditioned temperature signals and the signal 
conditioned load cell signals of ICt op , IC m id and ICbot, 
in addition to grounding and voltage feed lines. These 
conditioned output signals are current-driven instead of 
volt age- driven, leading to a superior noise suppression. 
From the slip ring, each electrical signal runs through 
a single shielded twisted-pair cable, again reducing the 
noise pick-up with respect to unshielded straight cables, 
and is acquired by data acquisition modules from Beck- 
hoff (Vcrl, Germany) operating at a maximum sampling 
rate of 1 kHz at 16-bit resolution. 

The shear stress CTA signals are fed to a liquid 
mercury- type slip ring from Mercotac (Carlsbad, USA). 
The liquid mercury inside the slip ring is used as a signal 
carrier and eliminates the electrical contact noise inher- 
ent to brush-type slip rings. This feature is important 
as additional (fluctuating) resistance between the CTA 
probe and the CTA controller can drastically reduce the 
accuracy of the measurement. 

The sensors embedded in the OC are also wired 
through the OC drive axle and exit to a custom-made 
slip ring from Moog (Boblingen, Germany). Those sen- 
sors include optional CTA probes and optical fiber sen- 
sors for two-phase flow. 

The T 3 C system is controlled with a combination 
of programmable logic controllers (PLC) from Beckhoff 
which interact with peripheral electronics, and with a PC 
running a graphical user interface built in National In- 



struments LAB VIEW that communicates with the PLCs. 
The controlled quantities include rotation of the cylin- 
ders, bubble injection rate, temperature, and inclination 
angle of the system. 



D. Rotation rate control 

The initial maximum rotation rates of the IC and OC 
are 20 Hz and 10 Hz, respectively. As long as the vibra- 
tion velocities occurring in the main ball bearings of the 
T 3 C fall below the safe threshold value of 2.8 mm/s rms, 
there is room to increase these maximum rotation rates 
in case the total torque on the drive shaft is still below 
the maximum torque of 200 Nm. The measured vibra- 
tion velocity in the T 3 C system was found to be much less 
than that threshold, even when the system was running 
at these initial maximum rotation rates. The system is 
capable of operating at even higher rotation rates. 

Both cylinders are driven by separate but identical ac 
motors, using a timing belt with toothed pulleys in the 
gear ratio 2:1 for the IC and 9:4 for the OC. Each motor, 
type 5RN160L04 from Rotor (Eibergen, Netherlands), is 
a four pole 15 kW squirrel-cage induction motor, pow- 
ered by a high-frequency inverter, Leroy Somer Unidrive 
SP22T, operating in closed-loop vector mode. A shaft en- 
coder mounted directly onto the ac motor shaft provides 
the feedback to the inverter. Electronically upstream of 
each inverter is a three-phase line filter and an electro- 
magnetic compatibility filter from Schaffner (Luterbach, 
Switzerland), types FN3400 and FS6008-62-07, respec- 
tively. Downstream, i.e. between the inverter and the ac 
motor, are two additional filters from Schaffner, FN5020- 
55-34 and FN5030-55-34, which modify the usual pulse- 
width-modulated driving signal into a sine-wave. Im- 
portant features of this filter technique are foremost the 
reduction of electromagnetic interference in the lab and 
the reduction of audible motor noise usually occurring in 
frequency-controlled motors. 

The rotation rate of each cylinder is independently 
measured by magnetic angular encoders, ERM200, from 
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FIG. 4: The measured instantaneous rotation rate fluctuations normalized by the time-averaged rotation rate. Here the inner 
cylinder was preset to rotate at 5 Hz. 



Heidenhain (Schaumburg, USA), mounted onto the di- 
rect drive shafts of the IC and the OC. The magnetic 
line count used on the IC is 1200, resulting in an angular 
resolution of 0.3°, and 2048 lines, resulting in an angular 
resolution of 0.18° on the OC. The given angle resolu- 
tions do not take into account the signal interpolation 
performed by the Heidenhain signal controller, improv- 
ing the resolution by a factor of 50 at maximum. Fig. 0] 
shows one measured time series of the rotation frequency 
when the system was rotating at (/ic) = 5 Hz. The 
figure shows the measured instantaneous rotation (/ic) 
rate fluctuation as a function of time. It is clearly shown 
that the rotation is stable within 0.01% of its averaged 
value ((/ic))- 



E. Temperature control 

The amount of power dissipated by degassed water at 
21 °C, in the case of a stationary OC and the IC rotat- 
ing at 20 Hz with smooth unaltered walls, is measured to 
be 10.0 kW. Without cooling this would heat the 0.111 
m 3 of water at a rate of « 1.3 K/min. As the viscosity 
of water lowers by 2.4% per K, it is important to keep 
the temperature stable (to within at least 0.1 K) to ex- 
clude viscosity fluctuations and thus errors. While the 
temperature- viscosity relation of water is well-tabulated 
and will be corrected for during measurements, other flu- 
ids like glycerin solutions might not share this feature. In 
the case of glycerin, a 1 K temperature increase can lower 
the viscosity by ~ 7%. 

To ensure a constant fluid temperature inside the TC 
tank, a 20 kW Neslab HX-750 chiller (Thermo Fisher 



Scientific Inc., Waltham, USA) with an air-cooled com- 
pressor and a listed temperature stability of 0.1 K is con- 
nected to the T 3 C. Two rotary unions, located at the 
bottom of the OC drive axle, are embedded inside each 
other and allow the coolant liquid to be passed from the 
stationary lab to the rotating OC, as shown in Fig. [5] 
Both the stainless steel bottom and top plate of the OC 
contain internal cooling channels which are covered by a 
5 mm-thick nickel-coated copper plate, which in turn is 
in direct contact with the TC's inner volume. 

The temperature inside the TC tank is monitored by 
three PT100 temperature sensors, each set up in a 4- 
leads configuration with pre-calibrated signal condition- 
ers IPAQ-H pius from Inor (Malmo, Sweden) with an ab- 
solute temperature accuracy of 0.1 K. The relative ac- 
curacy is better than 0.01 K. Each sensor is embedded 
at mid- height inside the wall of the hollow ICt op , IC m id, 
and ICbot sections, respectively (shown in Fig. [5]). Thus 
one can check for a possible axial temperature gradient 
across the IC. The sensors do not protrude through the 
wall so as to keep the outer surface smooth, leaving 1 
mm of stainless steel IC wall between the sensors and 
the fluid. Except for the direct contact area with the 
wall, the PTIOOs are otherwise thermally isolated. In- 
side each IC section, the signal conditioner is mounted 
and its electrical wiring is fed through the hollow drive 
axle, ending in an electrical slip ring. The average over 
all three temperature sensors is used as feedback for the 
Neslab chiller. An example of temperature time tracers 
is plotted in Fig. [6j which shows that the temperature 
stability is better than 0.1 K. This data was acquired 
with the IC rotating at 20 Hz, a stationary OC and wa- 
ter as the working fluid, resulting in a measured power 
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dissipation by the water of 10.0 kW. To check the effects 
of this small temperature difference on the TC flow, we 
calculate the Rayleigh number based on the temperature 
difference of A = 0.1 K over the distance {Lrb = 0.366 
m) of the middle and top sensor positions: the result 
is Ra = (3gLft B A/ ' nv — 5.9 xlO 6 (/? is the thermal ex- 
pansion coefficient, k the thermal diffusivity and v the 
kinematic viscosity). The corresponding Reynolds num- 
ber is estimated to be around Re^s ~ 0.25 x Ra ' 49 = 
500 [H, which is significantly smaller than the system 
Reynolds number of 2 x 10 6 . The effects of this small 
temperature gradient can thus be neglected in this high- 
Reynolds-number turbulent flow. 



F. Torque sensing 

Each of the ICbot, IC m id and ICt op sections are basi- 
cally hollow drums. Each drum is suspended on the IC's 
drive axle by two low-friction ball bearings, which are 
sealed by rubber oil seals pressed onto the outsides of 
the drums, encompassing the drive axle. A metal arm, 
consisting of two separate parts, is rigidly clamped onto 
the drive axle and runs to the inner wall of the IC section. 
The split in the arm is bridged by a parallelogram load 
cell (see Fig. [7]) . The load cells can be replaced by cells 
with different maximum-rated load capacity to increase 
the sensitivity to the expected torque. At this moment 
two different load cells, type LSM300 from Futek (Irvine, 
USA) , are in use with a maximum-rated load capacity of 
2224 N and 222.4 N, respectively. Each load cell comes 
with a pre-calibrated Futek FSH01449 signal conditioner 




FIG. 5: Left: Schematic sketch of the coolant flow through the 
T 3 C. Coolant enters at the bottom rotary union (blue) and 
flows straight up to the top plate by piping on the outside of 
the TC tank. Then the coolant enters the curling channels in 
the top plate (red) and is fed downwards again to run through 
the curling channels of the bottom plate before it exits at the 
second rotary union. Right: Bottom plate with the copper 
cover removed. 



18.10 




t (minutes) 



FIG. 6: Time traces of the measured temperature at three 
positions. The data were taken with the IC rotating at 20 
Hz, a stationary OC and water as the working fluid. The 
measured power dissipation by the water was 10.0 kW. 




FIG. 7: Horizontal cut-away showing the load-cell construc- 
tion inside the IC m id drum. The load cell spans the gap in 
the arm, connecting the IC drive axle to the IC wall. 



operating at 1 kHz, which is also mounted inside the 
drum. The electrical wiring is fed through the hollow 
drive axle to a slip ring on top. The hysteresis of each 
load cell assembly is less than 0.2 Nm, presented here 
as the torque equivalent. Calibration of the load cells is 
done by repeated measurements, in which a known se- 
ries of monotonically increasing or decreasing torques is 
applied to the IC surface. The IC is not taken out of 
the frame and is calibrated in situ. The torque is applied 
by strapping a belt around the IC and hanging known 
masses on the loose end of the belt, after having been 
redirected by a low-friction pulley to follow the direction 
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of gravity. 

Local fluctuations in the wall-shear stress can be mea- 
sured using the flush-mounted hot-film probes (type 
55R46 from Dantec Dynamics) on the surfaces of the in- 
ner and outer cylinder, shown in Fig. [5] 

An important construction detail determines how the 
torque is transferred to the load cell. Only the azimuthal 
component is of interest and the radial and axial compo- 
nents, due to possible non-azimuthal imbalances inside 
the drum or due to the centrifugal force, should be ig- 
nored. This is accomplished by utilizing the parallelo- 
gram geometry of the load cell, lying in the horizontal 
plane. It is evident that during a measurement the ro- 
tation rate of the IC should be held stable to prevent 




the rotational inertia of the IC sections from acting as a 
significant extra load. The mass balance of the system is 

important in decreasing this effect FIG. 8: A sketch of local sensors and bubble injectors in the 

system. 



G. Balancing and vibrations 

All three sections of the IC and the entire OC are sepa- 
rately balanced, following a one-plane dynamical balanc- 
ing procedure with the use of a Smart Balancer 2 from 
Schenck RoTec (Auburn Hills, USA). The associated ac- 
celerometer is placed on the main bottom ball bear- 
ing. The balancing procedure is reproducible to within 
5 grams leading to a net vibration velocity of below 2 
mm/s rms at the maximum rotation rates. According to 
the ISO standard 101816-1 [58[ regarding mechanical vi- 
brations, the T 3 C system falls into category I, for which 
a vibration velocity below 2.8 mm/s rms is considered 
acceptable. 

Another feature is the air springs, i.e. pressure regu- 
lated rubber balloons, placed between the floor and each 
of the four support feet of the T 3 C frame. They lift the 
frame fully off the ground and hence absorb vibrations 
leading to a lower vibration severity in the setup itself 
and reducing the vibrations passed on to the building. 

Two permanently installed velocity transducers from 
Sensonics (Hertfordshire, UK), type PZDC 56E00110, 
placed on the top and bottom ball bearings, constantly 
monitor the vibration severity. An automated safety 
PLC circuit will stop the IC and OC rotation when 
tripped. Thus, dangerous situations or expensive repairs 
can be avoided, as this acts as a warning of imminent 
ball bearing failure or loss of balance. 



H. Bubble injection and gas concentration 
measurement 

Eight bubble injectors, equally distributed around the 
outer perimeter of the TC gap as shown in Fig. [FJ are 
built into the bottom plate of the TC tank. Each bubble 
injector consists of a capillary housed inside a custom- 
made plug ending flush with the inside wall. They can 
be changed to capillaries of varying inner diameter: 0.05, 



0.12, 0.5 and 0.8 mm. This provides a way to indirectly 
control the injected bubble radius, estimated to be in 
the order of 0.5 mm to 5 mm, depending on the shear 
stress inside the TC. Smaller bubbles of radius less than 
0.5 mm, or microbubbles, can be injected by replacing 
the capillaries inside the plugs with cylinders of porous 
material. 

Two mass flow controllers from Bronkhorst (Ruurlo, 
Netherlands), series EL-Flow Select, are used in parallel 
for regulating the gas, i.e. filtered instrument air, with 
a flow rate at a pressure of 8 bars. One controller with 
a maximum of 36 1/min takes care of low-gas volume 
fractions, and the second controller with a maximum of 
180 1/min of high-gas volume fractions, presumably up 
to 10%. The gas enters the gap by a third rotary union 
located at the very bottom of the OC drive axle, below 
the coolant water rotary unions. Thus the OC drive axle 
has three embedded pipes running through its center, 
which are fed by three rotary unions at the bottom. Each 
pipe is split into separate channels again inside of the OC 
drive assembly to be routed where needed. 

Vertical channels running through the near-center of 
the TC tank's top plate connect the tank volume to a 
higher located vessel in contact with the ambient air. 
Excess liquid or gas can escape via this route to prevent 
the build-up of excessive pressure. We refer to it as an 
expansion vessel. The expansion vessel can also be used 
to determine the global gas volume fraction inside the 
TC tank. The vessel is suspended underneath a balancer 
which continuously registers the vessel's mass. Starting 
at zero percent gas volume fraction to tare the vessel's 
mass, one can calculate the global gas volume fraction 
by transforming the liquid's mass that is subsequently 
pushed into the vessel by the injected gas, into its equiv- 
alent volume. In the case of a rotating OC the stationary 
expansion vessel can not (yet) be connected and excess 
liquid is collected in a stationary collecting ring encom- 
passing the rotating top plate. 
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The original concept for measurement of the global gas 
volume fraction would not have had the restriction of 
a stationary OC. It makes use of a differential pressure 
transducer attached to the OC that measures the pres- 
sure difference between the top and bottom of the TC 
tank. Comparing this difference to the expected single- 
phase hydrostatic pressure difference, one could calculate 
the gas volume fraction. This method depends on the 
dynamic pressure being equal at the top and bottom. It 
fails however, due to unequal dynamic pressure induced 
by secondary flows. This is not unexpected as a wide 
variety of flow structures can exist in turbulent TC flow, 
like Taylor- vortices. 



I. Optical access and local sensors 

Flow structure and velocity fluctuations of TC flow 
have been studied extensively at low Reynolds num- 
bers, but few experiments have been performed at high 
Reynolds number (Re > 10 5 ). Previous velocity mea- 
surements were mainly done with intrusive measurement 
techniques like hot-film probes. Indeed it was found that 
the wake effects induced by an object inside a closed ro- 
tating system can be very strong (59[. Better velocity 
measurements inside the TC gap use nonintrusive opti- 
cal techniques such as LDA [38[, PIV [H, PTV. The 
optical properties of the outer cylinder are hence crucial 
for this purpose. 

The outer cylinder of the T 3 C is transparent, and four 
small areas of the top and bottom plates consist of view- 
ing portholes made of acrylic to allow for optical access 
in the axial direction, as can be seen in Figs. [5] and [7] 
The outer cylinder was thermally treated to homoge- 
nize the refractive index and to remove stresses inside 
the acrylic. Thanks to this optical accessibility, all three 
velocity components inside the gap can be measured opti- 
cally. For the velocity profile measurements, we use LDA 
(see Sec. HUH). 

Various experimental studies have been done to ex- 
amine bubbly DR. However, two main issues of bubbly 
DR in turbulent TC are still not well-studied. How do 
bubbles modify the turbulent flow? How do bubbles dis- 
tribute and move inside the gap? It is certainly impor- 
tant to measure local liquid and bubble information in- 
side the gap. Various local sensors, shown in Fig. [5J can 
be mounted to the T 3 C system. Here we highlight two of 
them: the phase-sensitive CTA and the 4-point optical 
fiber probe. 

1. Phase-sensitive CTA 

Optical techniques (such as LDA and PIV) are only ca- 
pable of measuring flow velocities in a bubbly flow when 
the gas volume fraction is very low (typically less than 
1%). Hot-film measurements in bubbly flows also im- 
pose considerable difficulty due to the fact that liquid 



and gas information is present in the signal. The chal- 
lenge is to distinguish and classify the signal correspond- 
ing to each phase. The hot-film probe does not provide 
by itself means for successful identification [(3(| ■ To over- 
come this problem, a device called phase-sensitive CTA 
has been developed (see Refs. |6ll - [64| ). In this technique, 
an optical fiber is attached close to the hot-film so that 
when a bubble impinges on the sensor it also interacts 
with the optical fiber. The principle behind the optical 
fiber is that light sent into the fiber leaves the fiber tip 
with low reflectivity when immersed in water, and with 
high reflectivity when immersed in air. Hence the fiber is 
able to disentangle the phase information by measuring 
the reflected light intensity. It has proved to be an use- 
ful tool for liquid velocity fluctuation measurements in 
bubbly flows [63[ . Phase-sensitive CTA probes are only 
mounted through the holes of the outer cylinder when 
necessary. 

2. J^-Point optical probe 

Instead of using a single optical fiber to discriminate 
between phases as described in the previous paragraph, 
one can construct a probe consisting of four such fibers. 
The four fiber tips are placed in a special geometry: three 
fiber tips of equal length are placed parallel in a trian- 
gle, and the fourth fiber is placed in the center of gravity 
and protrudes past the other fiber tips (see Ref. [65| for 
a schematic of the probe). Knowing this geometry and 
processing the four time series on the reflected light in- 
tensity, it becomes possible to estimate not only the size 
of the bubble that impinges onto the fiber tips, but also 
the velocity vector and the aspect ratio. To measure the 
bubble distribution inside the TC gap and other bubble 
dynamics, the 4-point optical probe is mounted through 
the holes of the outer cylinder only when necessary. We 
refer to Ref. (6(| for details on the measurement princi- 
ple of the 4-point optical probe. Support for this probe is 
built into the T 3 C by incorporating opto-electrical con- 
verters into the outer cylinder's bottom plate. 

Other probes used, such as temperature and shear 
stress sensors will not be described as they are standard. 

III. EXAMPLES OF RESULTS 

In this section we will demonstrate that the facility 
works by outlining our initial observations of the torque, 
velocity profiles and bubbly effects. 

A. Torque versus Reynolds number for single 
phase flow 

We first measure the global torque as a function of 
the Reynolds number in the present T 3 C apparatus with 
a stationary outer cylinder. The torque (T m ;d) on the 
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FIG. 9: (a) The non-dimensional torque and (b) the compen- 
sated torque GRe" 1 ' 80 versus Reynolds number in the high- 
Reynolds-number regime for the measurements with the T-M 
TC (open triangles) and T 3 C (open circles) apparatuses. 



FIG. 10: (a) The frequency regime of counter-rotation in the 
present measurements, (b) The dimensionless torque G versus 
inner cylinder Reynolds number Rei for the different outer 
cylinder frequencies. 



middle section of the inner cylinder is measured for 
Reynolds numbers varying from 3 x 10 5 to 2 x 10 6 . 
We use the same normalization as Ref. [3J) to define the 
non-dimensional torque as G = T m id/27T/Cw 2 L m id. The 
present measurements are performed in the T 3 C with a 
radius ratio r\ = 0.716, which is close to the value 77 = 
0.725 of the T-M TC examined by Lathrop et al. 
Figure^a) shows G versus Re^ for Re^ > 3x 10 5 . The fit- 
ting exponent for the data by Lathrop et al. [35j is 1.86, 
and the result of the present measurement gives 1.75. 

To better compare these data, we use G -180 to com- 
pensate for the torque GRe^ 1 ' 80 , which is shown in Fig. 
[9]Jb). Overall, the data of Ref. (HJ shows a higher ex- 
ponent than the compensated value of 1.80, and the 
present measurement is lower. However, both data sets 
clearly exhibit deviations from a single power law. The 
present measurement shows an oscillating trend when the 
Reynolds number is higher than 8 x 10 5 . This is likely 
induced by transitions between different flow structures, 
which will be studied systematically in the T 3 C appara- 
tus with high scrutiny as this trend was not anticipated. 

We also examine the torque versus Reynolds number 
in the counter-rotation regime. The parameter space for 



the present measurements is shown in Fig. UOf a). For 
three fixed outer cylinder frequencies f — -3, -6, and -10 
Hz, we measure the torque on the inner cylinder with 
increasing inner cylinder frequency. Fig. UOf b) shows 
the measured G versus Re^ for different outer cylinder 
frequencies. The results of pure inner rotation are also 
plotted in the figure for comparison. As shown in Fig. 
fTUT b). both the slope and the amplitude of G versus Re; 
for counter-rotation are different from those of pure inner 
rotation. It is not unexpected that the torque amplitude 
of the counter-rotation case exceeds the pure inner rota- 
tion case. The data for f„ < are not simple translations 
of the pure inner cylinder curve, but depend on the rota- 
tion frequency of the outer cylinder in a non-trivial way. 
In Ref. [56| we offer a unifying view to better under- 
stand G(Re;, Re ) in the regime of large inner and outer 
Reynolds numbers. 



B. Velocity profile measured with LDA 

Previous velocity measurements were mainly done with 
intrusive measurement techniques like hot-film probes. A 
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better way for velocity measurements inside the TC gap 
are nonintrusive optical techniques such as LDA. Due to 
the curvature, the angle between two LDA beams slightly 
depends on the radial position along the gap. We have 
corrected for this via a ray-tracing calculation based on 
the system parameters [67j. Refraction effects have also 
been taken into account. To test the reliability of the cor- 
rection, we perform a velocity profile measurement when 
the system is in solid-body rotation, i.e. / = fa = f a = 2 
Hz. Open circles in Fig. [Til show the azimuthal velocity 
profile measured with LDA. The exact solid-body veloc- 
ity profile, r/r a , is shown with the solid line in the figure. 
The velocity has been normalized by its value on the in- 
ner wall of the outer cylinder. Figure [TT] clearly indicates 
that the LDA measurements agree with the solid-body 
profile (within 0.6%). 

When only the outer cylinder rotates, the TC flow with 
an infinite aspect ratio is linearly stable. The laminar 
azimuthal Coucttc velocity profile for infinite aspect ratio 
(i.e. without end-plate effects) reads [68| : 



viam(r) = Ar + B/r,A 



27Tf 

l-r] 2 



,B 



1 -rj 2 



(1) 



The corresponding laminar profile with the present sys- 
tem parameters, for f Q — 7.2 Hz, is shown with the 
dashed line in Fig. 1111 The open squares are the mea- 
sured azimuthal velocity profile with LDA when the outer 
cylinder is rotating at f a = 7.2 Hz. The measurement 
is carried out with the T 3 C at a statistically stationary 
state. The velocity has been normalized by its value on 
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FIG. 11: Open circles: the azimuthal velocity profile mea- 
sured with LDA when the system is in solid-body rotation 
(/« — fo — 2 Hz). The solid line corresponds to r/r of the 
solid-body flow. Open squares: the azimuthal velocity pro- 
file measured with LDA in the T 3 C system with a stationary 
inner cylinder (fi = 0, f a = 7.2 Hz). The dashed line corre- 
sponds to the laminar profile described in Eqn. |T]), applicable 
to a laminar profile between the cylinders, with no top and 
bottom plate effects. 



the inner wall of the outer cylinder. A clear deviation is 
found between the measurement results with the laminar 
profile. The reason for this deviation is due to end-plate 
effects H3, [H, [lOl - the laminar flow profile of type (HJ 
only exists when the aspect ratio of TC is infinite. More 
measurements are upcoming for studying velocity pro- 
files in high-Reynolds-number Taylor- Couette flow. Cer- 
tainly the end-plate effects will be significantly reduced 
for the turbulent cases, when the liquid velocity fluctu- 
ations dominate the flow. Additional results on velocity 
profiles in T 3 C and a detailed discussion will be published 
elsewhere. 



Torque versus gas concentration for bubbly TC 
flow 



The T 3 C apparatus was designed for studying both 
single- and two-phase flows. We perform torque measure- 
ments for a bubbly flow in the present T 3 C apparatus 
with a stationary outer cylinder. In these first experi- 
ments with the T 3 C facility for bubbly flow we compare 
the torque with the previous measurements in T-M TC 
by van den Berg and coworkers [53, 54]. The air bubbles 
are injected into the turbulent flow through eight capil- 
laries located at the bottom plate. The gas fraction a is 
measured by weighing the excess volume pushed out of 
the system because of the added gas [H, [54| . The bub- 
ble radius is dictated by the strength of the shear and 
is typically in the range of 3—0.2 mm, which is kept as 
similar to Ref. [HJ . For three constant inner cylinder ro- 
tation frequencies, we measure the torque whilst slowly 
increasing the void fraction. 

When injecting bubbles into the flow, the effective 
kinematic viscosity [7l| and density are changed as fol- 
lows (for low gas concentrations): 



Pa = p(l 



(2) 
(3) 



where a is the gas volume concentration, and p and v are 
the density and viscosity for a = 0. The non-dimensional 
torque is defined based on the two-phase viscosity and 
density as: 



G(a) 



T 



mid 



2wp a ulL mid 271-/9(1 - a)v 2 (1 + § a) 2 L 



mid 



(4) 

with T m ;d as the measured torque on the middle section 
of the inner cylinder. We focus on the effect of bubbles 
on the torque and neglect the change of the Reynolds 
number. 

Figure [T2] shows the measured torque for two-phase 
flow compensated for by the single-phase torque as a 
function of the gas volume concentration. The linked 
dots are the present measurement results for different 
Reynolds numbers. As shown in the figure, a reduction 
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of gas-void fraction on non-dimensional tor que: dashed-lines 
are results measured in the T-M TC system [541 ]: linked-solid 
dots are measurement results in the present T 3 C system. 



in the required torque for constant rotation rate is ob- 
served with increasing void fraction. The drag reduction 
is generally larger than the density and viscosity effects 
on the torque, as shown with the solid line in the figure. 
The measured maximum drag reduction increases with 
increasing Reynolds number, which is consistent with 
Ref. [53| . As shown in the figure, more than 50% of drag 
reduction has been achieved by only adding 4% gas into 
the system for Re; = 2 x 10 6 , which was not attainable 
by previous measurements in other TC apparatuses. 

For quantitative comparison, the measurements in T- 
M TC [54], with Re = 5.5 x 10 5 and 7.9 x 10 5 , are shown 
as the dashed lines in the figure. The agreement is clearly 
revealed by the nice collapse of two data sets: the data 
of Re = 5.1 x 10 5 in T 3 C and the data of Re = 5.5 x 10 5 
in T-M TC. The data of Re = 7.9 x 10 5 in [H also lies 
between the present data of Re = 5.1 x 10 5 and 1.0 x 10 6 . 
Thus good agreement for bubbly drag reduction has been 
found between the T-M TC and T 3 C systems. 

Note that the new T 3 C system shows surprisingly large 
drag reduction (more than 50%) thanks to the system's 
capacity for high Reynolds numbers. This tremendous 
drag reduction in the high-Reynolds-number regime is 
very important to better understand the mechanism of 
turbulent drag reduction in bubbly flows, since large 
Reynolds numbers are more relevant to applications (e.g., 
carrier ships). More measurements will be done in this 
unexplored parameter space in the near future, and in 
particular we will explore the spatial bubble distribution. 
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IV. SUMMARY AND OUTLOOK 

A new turbulent Taylor-Couette apparatus, named 
T 3 C, has been developed, consisting of two independently 
rotating cylinders. The torque measurements for pure in- 
ner cylinder rotation agree well with previous results in 
the overlapped parameter regime. We also performed 
experiments in the unexplored parameter regime of large 
Reynolds number and counter-rotation of the cylinders. 
The data of G versus Re; for the counter-rotation situa- 
tion are not simple translations of the pure inner cylin- 
der curve, but depend on the rotation frequency of the 
outer cylinder in a nontrivial way. The nonintrusive 
measurement technique LDA is applied to the system 
for measuring velocity profiles through the transparent 
outer cylinder of the T 3 C, and it has proved to be an 
excellent tool for flow velocity measurements inside TC. 
The torque measurements of bubbly flow in the T 3 C sys- 
tem showed surprisingly large drag reduction (more than 
50%) at high Reynolds numbers, which was not attain- 
able by previous measurements in other TC apparatuses. 

The inner cylinder Reynolds number Re; of older ex- 
isting turbulent TC facilities can be as high as 10 6 . How- 
ever, another two control parameters, the outer cylinder 
Reynolds number Re D and radius ratio rj, still have not 
been explored in the highly turbulent regime (Re > 10 5 ). 
With this newly-built Taylor-Couette system, the param- 
eter space (Rej, Re Q , rf) has been extended to (2.0 x 10 6 , 
±1.4 x 10 6 , 0.716 - 0.909), as shown in Fig. EES Various 
research issues of single-phase turbulent TC flows can be 
studied with this new T 3 C facility; for example, turbu- 
lent momentum transport in the unexplored parameter 
space of co- and counter-rotation, the role of the "geo- 
metrical quasi-Prandtl number" |3[ on turbulent momen- 
tum transport, angular momentum/velocity profiles, and 
the connections with the global torque, and statistics of 
turbulent fluctuations of different velocity components. 
These studies will bridge the gap between RB flow and 
TC flow in the turbulent regime toward a better under- 
standing of closed turbulent systems. 
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FIG. 14: A snapshot of the bubble distribution in the T 3 C, 
Rei = 1.0 x 10 6 . The outer cylinder is stationary. 



The injection of bubbles offers another degree in the 
parameter space for studying two-phase TC flows. The 
study of the bubbly DR can be improved further by com- 
bining the global torque and the local measurements of 
the bubble velocity and distribution (see Fig. [T4|) . Apart 
from studying bubbly DR, T 3 C is also an ideal system for 



[1] P. Bradshaw, J. Fluid Mech. 36, 177 (1969). 
[2] B. Dubrulle and F. Hersant, Eur. Phys. J. B 26, 379 
(2002). 

[3] B. Eckhardt, S. Grossmann, and D. Lohse, J. Fluid Mech. 

581, 221 (2007). 
[4] G. Ahlers, S. Grossmann, and D. Lohse, Rev. Mod. Phys. 

81, 503 (2009). 
[5] D. Lohse and K. Q. Xia, Annu. Rev. Fluid Mech. 42, 335 

(2010). 

[6] M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65, 
851 (1993). 

[7] E. Bodenschatz, W. Pesch, and G. Ahlers, Annu. Rev. 

Fluid Mech. 32, 709 (2000). 
[8] B. Castaing, G. Gunaratne, F. Heslot, L. Kadanoff, A. 

Libchaber, S. Thomae, X. Z. Wu, S. Zaleski, and G. 

Zanetti, J. Fluid Mech. 204, 1 (1989). 
[9] T. Gebhardt and S. Grossmann, Z. Phys. B. Condens. 

Matter 90, 475 (1993). 
[10] X. Chavanne, F. Chilla, B. Castaing, B. Hebral, B. 

Chabaud and J. Chaussy, Phys. Rev. Lett. 79, 3648 

(1997). 

[11] J. A. Glazier, T. Segawa, A. Naert, and M. Sano, Nature 
398, 307 (1999). 

[12] J. Niemela, L. Skrbek, K. R. Sreenivasan, and R. Don- 
nelly, Nature 404, 837 (2000). 

[13] E. Brown, D. Funfschilling, A. Nikolaenko, and G. 
Ahlers, Phys. Fluids 17, 075108 (2005). 

[14] C. Sun, L.-Y. Ren, H. Song, and K.-Q. Xia, J. Fluid 



studying other issues in dispersed multiphase flows; for 
example, the Lagran g i an p roperties of particles/bubbles 
in turbulence [63|, |72|, LX3(1 , and liquid velocity fluctua- 
tions induced by dispersed bubbles/particles [Til. |75|. It 
was found that even in homogeneous and isotropic tur- 
bulence, particles, drops and bubbles are not distributed 
homogeneously, but cluster [76]. One example of bub- 
ble clustering in turbulent TC flow is shown in Fig. 1141 
Given the rich flow structures inside turbulent TC flow, it 
is of great interest to quantify particles/bubbles cluster- 
ing with varying turbulent structures, and this will also 
be done in future research. 



Acknowledgments 

This study was financially supported by the Tech- 
nology Foundation STW of the Netherlands. The au- 
thors would like to thank S.G. Huisman for velocity pro- 
file measurements, G. Ahlers, R. Delfos, S. Grossmann, 
T. Mullin, D. Narezo, G.P. Pfister, E. van Rietbergen, 
R.J. A.M. Stevens, T. van Terwisga, and J. Westerweel 
for scientific discussions. We also gratefully acknowledge 
technical contributions by B. Benschop, S.J. Boorsma, 
M. Bos, G. Mentink, R. Nauta, R.O. Heuvel, S. Bek- 
mann, and J. Schepers. D.P. Lathrop was supported by 
a grant from the NSF/DMR of the USA. 



Mech. 542, 165 (2005). 

[15] R. du Puits, C. Resagk, A. Tilgner, F. H. Busse, and A. 
Thess, J. Fluid Mech. 572, 231 (2007). 

[16] G. Ahlers, E. Bodenschatz, D. Funfschilling, and J. Hogg, 
J. Fluid Mech. 641, 157 (2009). 

[17] R. J. A. M. Stevens, R. Verzicco, and D. Lohse, J. Fluid 
Mech. 643, 495 (2010). 

[18] See the Focus issue of New Journal of Physics on New 
Perspectives in High-Rayleigh-Number Turbulent Con- 
vection edited by J. Schumacher and E. Bodenschatz, 
January (2011). 

[19] R. C. DiPrima and H. L. Swinney, in Hydrodynamic In- 
stabilities and the Transition to Turbulence, edited by H. 
L. Swinney and J. P. Gollub (Springer, New York, 1981), 
p. 139180. 

[20] C. D. Andereck, S. S. Liu, and H. L. Swinney, J. Fluid 

Mech. 164, 155 (1986). 
[21] M. A. Dominguez-Lerma, D. S. Cannell, and G. Ahlers, 

Phys. Rev. A 34, 4956 (1986). 
[22] D. Coles, J Fluid Mech 21, 385 (1965). 
[23] G. I. Taylor, Philos. Trans. R. Soc. London, Ser. A 223, 

289 (1923). 

[24] G. Pfister and I. Rehberg, Phys. Lett. 83, 19 (1981). 
[25] T. Mullin, K. A. Cliffe, and G. Pfister, Phys. Rev. Lett. 

58, 2212 (1987). 
[26] G. Pfister, H. Schmidt, K. A. Cliffe, and T. Mullin, J. 

Fluid Mech. 191, 1 (1988). 
[27] P. Buchel, M. Luecke, D. Roth, and R. Schmitz, Phys. 



13 



13 



Published in Rev. Sci. Instrum. 82, 025105 (2011). 



Rev. E 53, 4764 (1996). [53 
A. Esser and S. Grossmann, Phys. Fluids 8, 1814 (1996). 
S. Dong, J Fluid Mech 587, 1 (2007). [54 
S. Dong, J. Fluid Mech. 615, 371 (2008). 
G. P. Smith and A. A. Townsend, J. Fluid Mech. 123, [55 
187 (1982). 

A. A. Townsend, J. Fluid Mech. 144, 329 (1984). [56 
F. Wendt, Arch. Appl. Mech. 4, 577 (1933). 

F. Ravelet, R. Delfos, and J. Westerweel, Phys Fluids 22, [57 
055103 (2010). 

D. P. Lathrop, J. Fineberg, and H. L. Swinney, Phys. [58 
Rev. Lett. 68, 1515 (1992). 

D. P. Lathrop, Ph.D. thesis, University of Texas at 
Austin, 1992. 

D. P. Lathrop, J. Fineberg, and H. L. Swinney, Phys. 
Rev. A 46, 6390 (1992). [59 

G. S. Lewis, Ph.D. thesis, University of Texas at Austin, 

1996. [60 
G. S. Lewis and H. L. Swinney, Phys. Rev. E 59, 5457 

(1999) . [61 
T. H. van den Berg, C. Doering, D. Lohse, and D. Lath- 
rop, Phys. Rev. E 68, 036307 (2003). [62 
S. Grossmann and D. Lohse, J. Fluid. Mech. 407, 27 

(2000) . [63 
S. Grossmann and D. Lohse, Phys. Rev. Lett. 86, 3316 

(2001) . [64 
S. Grossmann and D. Lohse, Phys. Rev. E 66, 016305 

(2002) . [65 
S. Grossmann and D. Lohse, Phys. Fluids 16, 4462 
(2004). [66 
F. H. Busse, Lect. Notes Phys. 12, 103 (1972). 

S. L. Ceccio, Annu. Rev. Fluid Mech. 42, 183 (2010). [67 
N. K. Madavan, S. Deutsch, and C. L. Merkle, Phys. [68 
Fluids 27, 356 (1984). 

T. Takahashi, A. Kakugawa, M. Makino, T. Yanagihara, [69 
and Y. Kodama, in Proceedings of the 74th General [70 
meeting of SRI, 2000. 

Y. Kodama, A. Kakugawa, S. Nagaya, and T. Kawamura, 
Report of 24th US- Japan Joint Meeting, Marine Facilities [71 
Panel of the U.S. -Japan Cooperative Program in Natural [72 
Resources, 2001. [73 
Y. Kodama, A. Kakugawa, T. Takahashi, and H. 
Kawashima, Int. J. Heat Fluid Flow 21, 582 (2000). [74 
C. Merkle and S. Deutsch, Frontiers in Experimental [75 
Fluid Mechanics, Lecture notes in Engineering Vol. 46, 
p. 291, edited by M. G. el Hak (Springer, Berlin, 1989). [76] 
R. Latorre, A. Miller, and R. Philips, Ocean Eng. 30, 
2297 (2003). 



T. H. van den Berg, S. Luther, D. P. Lathrop, and D. 

Lohse, Phys. Rev. Lett. 94, 044501 (2005). 

T. H. van den Berg, D. P. M. van Gils, D. P. Lathrop, 

and D. Lohse, Phys. Rev. Lett. 98, 084501 (2007). 

Y. Murai, H. Oiwa, and Y. Takeda, Phys. Fluids 20, 

034101 (2008). 

D. P. M. van Gils, S. G. Huisman, G. -W. Bruggert, C. 
Sun, and D. Lohse, Phys. Rev. Lett. 106, 024502 (2011). 
M. S. Paoletti, and D. P. Lathrop, Phys. Rev. Lett. 106, 
024501 (2011). 

Part 1: General Guidelines: Mechanical Vibration - 
Evaluation of Machine Vibration by Measurements on 
Non-rotating Parts (International Organization for Stan- 
dardization, Geneva, Switzerland, 1995), ISO 10816- 
1:1995. 

C. Sun, T. Mullin, L. van Wijngaarden, and D. Lohse, J. 
Fluid Mech. 664, 150-173 (2010). 

R. Zenit, D. L. Koch, and A. S. Sangani, J. Fluid Mech. 
429, 307 (2001). 

S. Luther, J. Rensen, T. H. van den Berg, and D. Lohse, 
Exp. Therm. Fluid Sci. 29, 821 (2005). 
T. H. van den Berg, Ph.D. dissertation, University of 
Twente, 2006. 

J. Martinez, D. Chehata Gomez, D. P. M. van Gils, C. 

Sun, and D. Lohse, J. Fluid Mech. 650, 287 (2010). 

T. H. van den Berg, W. D. Wormgoor, S. Luther, and D. 

Lohse, Macromol. Mater. Eng. (in press). 

S. Guct, R. V. Fortunati, R. F. Mudde, and G. Ooms, 

Part. Part. Syst. Charact. 20, 219 (2003). 

J. Xue, M. Al-Dahhan, M. P. Dudukovic, and R. F. 

Mudde, Flow Meas. Instrum. 19, 293 (2008). 

S. G. Huisman, M. S. thesis, University of Twente, 2010. 

L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perg- 

amon Press, Oxford, 1987). 

D. Coles and C. Van Atta, J. Fluid Mech. 25, 513 (1966). 
B. Dubrulle, O. Dauchot, F. Daviaud, P.-Y. Longaretti, 
D. Richard, and J.-P. Zahn, Phys. Fluids 17, 095103 
(2005). 

A. Einstein, Ann. Phys. 19, 289 (1906). 

D. Lohse, Physics 1, 18 (2008). 

F. Toschi and E. Bodenschatz, Annu. Rev. Fluid Mech. 
41, 375 (2009). 

M. Lance and J. Bataille, J. Fluid Mech. 222, 95 (1991). 
J. M. Rensen, S. Luther, and D. Lohse, J. Fluid Mech. 
538, 153 (2005). 

E. Calzavarini, M. Kerscher, D. Lohse, and F. Toschi, J. 
Fluid Mech. 607, 13 (2008). 



14 



14 



